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Abstract 

We describe the public ESO near-IR variability survey (VVV) scanning 
the Milky Way bulge and an adjacent section of the mid-plane where star 
formation activity is high. The survey will take 1929 hours of observations 
with the 4-metre VISTA telescope during five years (2010 - 2014), covering 
~ 10 9 point sources across an area of 520 deg 2 , including 33 known globular 
clusters and ~350 open clusters. The final product will be a deep near-IR at- 
las in five passbands (0.9-2.5 /mi) and a catalogue of more than 10 6 variable 
point sources. Unlike single-epoch surveys that, in most cases, only produce 
2-D maps, the VVV variable star survey will enable the construction of a 3-D 
map of the surveyed region using well-understood distance indicators such 



1 



as RR Lyrae stars, and Cepheids. It will yield important information on the 
ages of the populations. The observations will be combined with data from 
MACHO, OGLE, EROS, VST, Spitzer, HST, Chandra, INTEGRAL, WISE, 
Fermi LAT, XMM-Newton, GAIA and ALMA for a complete understanding 
of the variable sources in the inner Milky Way. This public survey will pro- 
vide data available to the whole community and therefore will enable further 
studies of the history of the Milky Way, its globular cluster evolution, and 
the population census of the Galactic Bulge and center, as well as the inves- 
tigations of the star forming regions in the disk. The combined variable star 
catalogues will have important implications for theoretical investigations of 
pulsation properties of stars. 

Key words: Surveys, Stars: variables: general, Galaxy: bulge, Galaxy: disk 
PACS: 95.80. +p, 97.30. -b, 98.35.Jk, 98.35.Ln 

1 Introduction 

The bulk of the stars, gas and dust in the Milky Way are confined to its bulge 
and plane. As a result, in these directions, the extinction and crowding are high, 
making any study of the inner structure of the Galaxy difficult. Knowing how the 
stellar populations are distributed within the Galaxy is essential for such studies 
and hence the main goal of the described survey. Traditional distance indicators 
have been used with various success in the past. The approach was to concen- 
trate on clear "windows", where optical surveys can be carried out (e.g., MACHO, 
OGLE, EROS). In this paper we describe the VISTA Variables in the Via Lactea 
(VVV) surveyLl, an ESO (European Southern Observatory) public near-IR variabil- 
ity survey. Its area includes the Milky Way bulge and an adjacent section of the 
mid-plane where star-formation activity is high. This survey will be conducted in 
the period 2010 - 2014 and will map the whole bulge systematically for multiple 
epochs. 

We plan to cover a 520 deg 2 area (Fig.Q} containing ~10 9 point sources. Our 
survey will give the most complete catalogue of variable objects in the bulge, with 
more than ~10 6 variables. Chief among them are the RR Lyrae, which are ac- 
curate primary distance indicators, and well understood regarding their chemical, 
pulsational and evolutionary properties. For the sake of space and coherence we 
concentrate on the RR Lyrae and the star clusters, noting that similar studies can 
be done for many of the other populations of variable objects. 

Earlier single-epoch near-IR surveys (e.g., COBE, 2MASS, GL IMPSE) have 
proven that the Galactic bulge is triaxial and boxy, and contains a bar (|Dwek et al. , 

'Detailed information about the VVV survey can be found in http://vvvsurvey.org/ 
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Figure 1: 2MASS map of the inner Milky Way showing the VVV bulge (solid 
box, -10° < / < +10° and -10° < b < +5°) and plane survey areas (dotted box, 
-65° < I < -10° and -2° < b < +2°). 



1995k lLopez-Corredoira et all |2005|; benjamin et all 120051) . Presently, the only 



model we have for the formation of boxy/barred bulges is through secular evo- 
lution of a pre-existing disk. This scenario is believed to be the dominant chan- 
nel of formation of bulges in late-type spirals (Sbc), whereas early-type spiral 
bulges (SO/Sa) show structural and kinematic evidence for an early, rapid col- 
lapse, which seems to be confi rmed by the old age of their stellar populations (e.g., 



Kormendy & Kennicutt, 2004). 



However, the best-studied spiral bulge, that of the Milky Way, is precisely the 
most problematic one to understand in this context. While its surf ace brightness 
shows a barred structure, i ts stellar population is predominantely old dKuijken & Rich , 



20021 : IZoccali et all I2003h and has a-element enhancement, characteristic of rapid 
formation. Nevertheless, the high mean age of the Bulge still leaves space for a 
small fracti on of young stellar objects (YSO) which have b een found in the inner 
Bulge (e.g.. ISchuller et all 120061 : lYusef-Zadeh et all l2009h . This is in agreement 



with the results of 



IZoccali et all d2006h which indicate that the chemical compo- 
sition of the bulge stars is different from that of both thin and thick-disk stars. 
Thus, the predictions from the formation of the Milky Way bulge through secular 
evolution of the disk seem to be in conflict w ith some key properties of its stellar 
population. Ho wever, [M elendez et alj (120081) recently published results that are in 
contradiction to IZoccali et al.1 (|2006) and show that bulge and disk stars are indis- 



3 



tinguishable in their chemical composition. Given that the near-IR colours depend 
strongly on metallicity, the VVV survey will help us to investigate the metallic- 
ity distribution in the survey region. Spectroscopic data (e.g., future APOGEE; 



Majewski et all. 120070 will provide additional a-element abundances. 



Our survey of the RR Lyrae in the Ga lactic bulge will allow us to map its 3-D 
structure (as shown by lCarney et all 1 19951) and will provide key information on the 
age o f its population , given that RR Lyrae stars are tracers of the old population 
(e.g., Catelanir2004b . 20091. and references therein). This will enable us to combine 
the ages of the stellar populations with their spatial distributions. We note that most 
single-epoch surveys only provide 2-D maps. With the present survey, the peak and 
width of the RR Lyrae distribution is expected to be measured with an accuracy of 
better than 0.01 mag, which is the required precision to determine the 3-D structure 
not only of the bulge, but also of the Sagittarius dw arf spheroidal galaxy (Sgr dSph) 
located behind the Milky Way (e.g., lAlardl l 1996) and included in our survey. 

At the same time, a comparison between the RR Lyrae (and type II Cepheids) 
in the field and in glo bular cluster s may hold precious information about the forma- 
tion of the bulge (e.g. Feast et all 120081) . Modern ACDM cosmology predicts that 
large galaxies such as the Milky Way formed by accretion of hundreds of smaller 
"protogalactic fragments", perh aps not unlike the progenitors of the present-day 
dwarf spheroidal satellites (e.g., Abadi et al. . 2003() . Interestingly, two very mas- 
sive globular clusters in the Galactic bulge, NGC 6388 and NGC 6441, have re- 
cently been suggested to be the rem nants of dwar f galaxies that were accreted in 
the course of the Galaxys history (IRee et all 120021) . These clusters might prove 
similar to the cases of M54 (NGC 6715), in the center of the Sgr dSph, which is 
currently being cannibalised by the Milky Way (llbata et all 1 1995b . and of a) Cen 
(NGC 5 139) , which has long been suspected to be the remnant nucleus of a dwarf 
galaxy (e.g.. lAltmann et all 120051 and references therein). Our proposed search 
for RR Lyrae and type II Cepheids in the Galactic bulge will reveal the presence 
of debris related to the accretion events that might have left behind NGC 6441 as 
remnant object. The latter is part of our survey. 

In order to understand the Milky Way's populations globally, it is necessary 
to survey the inner Galactic plane as well. Therefore, we will survey an adjacent 
region of the mid-plane and provide a Legacy Database and 3-D atlas of a large 
Population I (i.e. young and luminous stars) region. We have selected the re- 
gion -65° < I < -10° and \b\ < 2° (see Fig.Q]), where star-formation activity is 
high and for which there will be complementary optical, mid-IR, and far-IR data 
from VPHAS+, the Spitzer, GLIMPSE and MIPSGAL surveys, and from the all- 
sky AKARI and WISE survey. The addition of this region will also permit us to 
discriminate between various models of the inner Galactic structure which, besides 



the triaxial bulge, contain a long bar and a ring (e.g.. lL6pez-Corredoira et all 2007 
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and references therein) or not (e.g.. iMerrifieldl. 12004, and references therein). In- 
deed, the selected region includes the putative negative-longitude tip of the long 
bar (at / « -14°, \b\ < 1°), which has not yet been observed. 

The large survey area will allow several remaining astrophysical problems to 
be addressed. For example, the effect of the environment on star formation and in 
particular the initial mass function (IMF) at low masses is presently poorly known. 
This issue will be addressed statistically by observing hundreds of star-forming 
regions and cross-correlating the shapes of their luminosity functions with cluster 
density, the presence of high-mass stars, and galactocentric distance. For compar- 
ison, VVV survey will reach 1 mag deeper than UKIDSS Galactic Plane Survey 
(GPS), which overlaps with VVV in the region of -2° < I < +10°, \b\ < 2°. 
Other important parameters, such as velocity dispersion and metallicity, will be 
determined by spectroscopic follow-up observations. In addition, the luminosity 
function of the clusters themselves will be measured, for both star-forming clusters 
and more evolved open clusters. 

These issues cannot be addressed with optical surveys, owing to the high ex- 
tinction in the plane. The Spitzer data will be invaluable for detecting the most 
obscured high-mass protostars within star-forming regions. A near-IR survey will 
be more sensitive to all but the reddest objects, and the superior spatial resolution 
in these wavebands will be essential for resolving distant clusters and the crowded 
field populations. 



2 Technical description 

2.1 Telescope and instrument design 

The Visible and Infrared Survey Telescope for Astronomy (VISTA) is a 4m-class 
"wide-field" telescope located at ESO's Cerro Paranal Observatory in Chile, de- 
signed to conduct large-scale surveys of the southern sky at near-IR wavelengths 
(0.9 to 2.5 /im). The telescope has an altitude-azimuth mount, and quasi Ritchey- 
Chretien optics. An f/l primary mirror was designed together with Cassegrain- 
focus instrumentation to offer the best solution to the difficult problem of combin- 
ing a wide-field with good i mage quality, and resul ts in a physically large focal 



plane with an f/3.25 focus dMcPherson et al. , 2006 ). VISTA's active optics uses 



two low-order curvature sensors, which operate concurrently with science expo- 
sures, and a high-order curvature sensor. 

The telescope is equipped with a near-IR camera containing 67 million pixels 
(an array of 16 x 2048 x 2048 Raytheon VIRGO IR detectors) of mean size ft' 34 
and available broad-band filters atZYJHK s and a narrow -band filter at 1.18 fim. 
Given VISTA's nominal pixel size, the diameter of the field of view is 1.65 deg. The 
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| | | | VLT / HAWK-I 

Figure 2: Diagram showing the array of the sixteen detectors on the VISTA camera 
and the axis orientation used to shift the camera in each exposure to obtain the tiles. 
For comparison we show the crescent Moon over the VISTA camera and the fields 
of view of UKIRT/WFCAM, HST/NICMOS, VLT/ISAAC, and VLT/HAWK-I. 

point-spread function (PSF) of the telescope+camera system (including pixels) is 
designed to have a full width at half maximum (FWHM) of Off 51, not including the 
contribution of atmospheric turbulence. Seeing, and other weather-related statistics 
for Cerro Paranal, are given at ESO's "Astroclimatology of Paranal" web page^|. 
The VISTA site is expected to have similar conditions, which are well suited to our 
survey requirements. 

The 16 detectors in the camera are not buttable and are arranged as shown in 
Figured Each individual exposure produces a sparsely sampled image of the sky 
known as a "pawprint", covering an area of 0.599 deg 2 . 

To 'fill in' the gaps between the detectors to produce a single filled "tile" with 
reasonably uniform sky coverage, the minimum number of pointed observations 

2 http://www.eso.org/gen-fac/pubs/astclim/paranal/ 
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(with fixed offsets) required is six (three offsets in Y and two offsets in X). After 
six steps an area of 1.501 deg 2 on the sky, corresponding to one tile, is (almost) 
uniformly covered. 



2.2 Data reduction 

We will use the enhanced VISTA Data F low System^ (VDFS : lEmerson et allbooj ; 



Irwin et all 120041 : lHambly et all 120041) . It includes all basic data reduction steps: 



(i) removing instrumental signature (bias and dark frames, twilight, and dome 
flatfields, linearity, bad pixel maps, cross-talk, gain calibrations), merging paw- 
prints into tiles and calibrating photometrically and astrometrically; 

(ii) extracting source catalogues on a tile-by-tile basis; 

(iii) constructing survey-level products - stacked pixel mosaics, difference im- 
ages, and merged catalogues; 

(iv) providing the team with both data access and methods for querying and 
analyzing the data; and 

(v) producing virtual observatory (VO)-compliant data products for delivery to 
the ESO archive. 

Figure|3] shows a flow chart of the data processing. The pipeline products are: 
astrometrically corrected and photometrically calibrated tiles in each filte r used, 
confidence maps, and homogeneous object catalogues ( Cross et al. , 20091) . The 
pipeline records the processing history and calibration information of each file, 
including calibration files and quality control parameters. The Cambridge Astron- 
omy Survey Unit (CASU) component of the VDFS will be responsible for the basic 
pipeline processing and the first calibration, all done on a daily basis. 



2.3 Combination/image subtraction (archive) 

The "second"-order data processing requires access to larger sets of data to produce 
survey products. It is carried out by the Wide Field Astronomy Unit's (WFAU) 
VISTA Science Archive (VSA) in Edinburgh. The Science Archive contains only 
calibrated data and catalogues, and no raw data. The Science Archive is responsible 
for: 

(i) image stacking to produce combined and differenced tiles and source merg- 
ing; 

(ii) quality control: assessment of the data quality and filtering of the data that 
do not meet the established criteria for photometric and astrometric accuracy; 

3 The VISTA Data Flow System (VDFS) is a collaboration between the UK Wide Field Astron- 
omy Unit at Edinburgh (WFAU) and Cambridge Astronomy Survey Unit (CASU), coordinated by 
the VISTA PI and funded for VISTA by the Science and Technology Facilities Council. 
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Figure 3: Flow chart of the VVV data processing (QC: Quality Control; ZP: Zero- 
point). 



(iii) light-cu rve extraction: this will be don e by im plementing an image-subtrac- 
tion algorithm ( Alard & Luptorl 19981 : Alard . 2000h . which will allow us to cre- 
ate the catalogue of variable sources. This method provides excellent results for 
crow ded fields in which the traditional aperture or PSF-fitting photometry fails 
(e.g., Ka luznv et all 120041 : Ipietrukowicz et al.ll2005h . 



2.4 Photometric calibration 

During the first period we will carry out the external calibrations and transforma- 
tions to the standard system using 2MASS and UKIDSS for bootstrapping^. 
The calibration scheme for a given filter is as follow: 

m ca i = m inst + ZP- k(X - 1) = m st _ + clr std (1) 

where m ca i is the calibrated magnitude, mi nst the measured instrumental magnitude, 
ZP the zeropoint, k the extinction coefficient, and X the airmass of the object. On 

4 The filter transmission curves for each instrument can be found at 
http://www.vista.ac.uk/index.html (VISTA), http://web.ipac.caltech.edu/staff/waw/2mass/opt_cal/in- 
dex.html (2MASS) and http://www.ukidss.org/technical/instrument/filters.html (UKIDSS). 
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the right-hand side of this equation, m st( i and c/r st d are the corresponding standard 
magnitude and colour. 

Calibration and quality control is done using 2MASS stars in the frames them- 
selves, applying col our equations to convert 2MASS photom etry to the VISTA 
photometric system ( Skrutskie et al. , 2006 : Hodgkin et al. , 20091) . 

There are thousands of unsaturated 2MASS stars in J H K s with photometric 
errors <0. 1 mag in every VISTA tile field. A large fraction of these can be suffi- 
ciently isolated even in the crowded fields. 

To calibrate the Y- and Z-band data (both filters are not available from 2MASS) 
we will use observations of the standard VISTA calibration fields as required by 
the ESO Public Survey Panel. Details will be published in a forthcoming paper 
describing the science verification. 

The internal gain correction applied through flat-fielding will place the detec- 
tors on a common zero-point system. After deriving this ZP in each tile, a double 
check using the overlap regions will be made to estimate the internal photometric 
accuracy. 



3 Observing strategy 



The VISTA tile field of view is 1.501 deg 2 , hence 196 tiles are needed to map the 
bulge area and 152 tiles for the dislH. Adding some X and Y overlap between tiles 
for a smooth match, the area of our unit tile covered twice is 1.458 deg 2 . Figured 
provides a schematic representation of the tiling scheme for the Galactic center 
region. 

The variability study in the bulge will be carried out in the K s band down to 
~18 mag (signal-to-noise « 3). The total exposure time for a VISTA tile field is 
162 s. Our strategy yields about 30 deg 2 per hour, or 300 deg 2 per night. The 
combined epochs will reach K s = 20 mag, which is three magnitudes fainter than 
the unreddened bulge main-sequence turn-off (MS turn-off), although the densest 
fields will be confusion-limited. However, applying both PSF fitting and image 
subtraction, we will recover the light curves of most objects down to K s = 18 mag, 
even in moderately crowded fields. This is more than 3 mag fainter than the unred- 
dened known RR Lyrae in the Galactic bulge. We expect to find RR Lyrae even in 
fields with Ay = 10 mag. 

Tabled] lists some reference ^ s -band magnitudes at the distance of the bulge 
for a range o f extinction and reddening value s. These typica l magn itud es were ob- 
tained from lCarney et a l. ( 1995), lAlarddl996l) . lAlcock et al1dl998l) . and lZoccali etal 



5 The tiles' spacing and orientation were calculated with the Survey Area Definition Tool (SADT) 
software to maximize the efficiency of the sky coverage. See http://www.vista.ac.uk/observing/sadt/ 
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Figure 4: Coverage of the Galactic center region overlaid on a mid-IR map. Boxes 
mark the tiles needed to cover the bulge area (see Fig.Q}, whereas the crosses mark 
the box centers. 



(2003). As a reference point, for Baade's window E(B - V) = 0.5 mag, so that 
Ay = 1.5, Aj = 0.4, and A K = 0.2 mag feieke & Lebofskvl. fl985h . This table 
shows that the tip of the bulge red-giant branch (RGB) will saturate (K s < 9.5), 
but for the RGB clump giants, and even for the tip of the RGB of the Sgr dSph 
galaxy, the VVV survey will be able to see giants throughout the bulge, even in 
the most obscured regions. The bulge RR Lyrae and the Sgr dSph galaxy red- 
clump giants will also be detected, even for the regions with the highest extinction 
(Ay > 30 mag) at low Galactic latitudes. Finally, the RR Lyrae of the Sgr dSph 
galaxy and the bulge MS turn-off stars will be detected only in the regions with 
low absorption (Ay < 10 mag) at higher latitudes. We are aware that the use of an 
'universal' extinction la w Ay=3.1 E(B — V) i s proble matic in the inner region of 
the bulge. According to iNishiyama et all (120061) and iGosling et all (12009b a single 
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extinction law is not consistent with the observations of the Galactic center along 
different lines of sight. 

For the plane survey, the Ag-band observations require a total time of 80 s on 
target, and an elapsed time of 366 s per tile. 

Bright point sources with K s < 9.5 mag will be saturated in the individual im- 
ages. This will therefore include most unreddened bulge Mira variables, but Miras 
in the Sgr dSph galaxy can be monitored, as well as Miras located in regions with 
very high extinction (e.g., next to th e Galactic center) . The Mira population in the 
Galactic center has been studied by Ma tsunag a et all d2009h . In addition, bright- 
star saturation may be an issue, but we estimate that even in the worst cases only a 
small portion of the field will be rendered useless. Hence we do not expect that the 
saturation of the brightest stars to effect our conclusions about three-dimensional 
structure of the inner Milky Way. For example, in the optical microlensing surveys 
where CCD bleeding is comparatively worse, less than 5% of the most crowded 
bulge fields are lost. 

To illustrate the precision of crowded-field IR photomeuy we include Fig.[5l 
showing photometry of the planetary transit OGLE-TR-113 obtained with NTT/ 
SOFI (top panel), and photometric accuracy of those observations as a function 
of magnitude (bottom panel). In order to evaluate the amplitude threshold for our 
detections, we have carried out Monte Carl o simulations using the R R Lyrae light 
curve templates from lJones et al. and lDel Principe et all (120051) . As a result, 

we find that, at a typical magnitude of K s « 15 - 16, and taking into account 
the expected photometric errors, we should be able to detect RR Lyrae stars with 
amplitudes down to Ak = 0.05 - 0.07 mag using 80 datapoints from the first three 
years of VVV operation, and further down to Ak = 0.03 - 0.05 mag if the dataset 
is extended to cover 1 80 phase points over a time frame of 5 years. 

The total estimated time per observing period is given in Table[2j which also 
shows the requirements for Moon, seeing, and transparency conditions. The times 
include overheads (both for readout and for changing to a new tile) and time pos- 
sibly spent on standard stars for the Z and Y observations (not used in the 2MASS 
and UKIDSS survey). This strategy allows us to provide various data to the com- 
munity, enabling follow-up throughout the survey. The full survey will require a 
total of 192 nights of observations over 5 years. A schematic schedule of the survey 
is shown in Fig. [6] 

During the first year, the whole bulge area will be observed in the ^ s band for 
6 consecutive epochs, for a total of 65 hours. A further 86 hours will be devoted 
to complete imaging of each bulge tile in Z Y J H. This will provide reliable near- 
simultaneous fluxes and colours for each tile area. 

The same strategy will be applied to the 152 tiles covering the disk area for 
the single-epoch and the quasi-simultaneous multi-colour disk survey. The total 
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Table 1: ^ ? -band magnitudes at the d istance of the bulge. The absorption values are based on the standard extinction law 
derived by iRieke & Lebofskvl dl985h . 
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Figure 5: Top panel: Light curve of the planetary transit OGLE-TR-113 measured 
with SOFI at the 3.5-m ESO NTT telescope in the ^ s band. This star is located 
in the Galactic plane in the field of Carina, a moderately crowded region. Bottom 
panel: Precision of the relative photometry obtained with SOFI as a function of 
magnitude. The magnitudes and rms are calculated in two iterations, removing the 
> 10cr outliers. 

time spent on the disk for the first year is thus 141 hours. Added to the 151 hours 
for the bulge, we will thus spend 292 hours in total on the survey during the first 
year. The multi-color observations, in combination with datasets from UKIDSS 
GPS (near IR), VST/VPHAS+ (optical) and GLIMPSE and GLIMPSE-II (mid- 
IR), will be used to build improved extinction maps for the survey region. Note 
that the individual, single-epoch observation blocks (OBs) in K s all have the same 
limiting depth (under the same conditions), whether forming part of the VVV's 
bulge or disk components. 

Being fully aware of the confusion and background limits, the observing plan 
would cycle alternately through fields of varying density for optimal sky subtrac- 
tion. The filter order in the OBs will be optimized to minimize overheads. 

During the second year, we will acquire another 20 epochs in K s for the whole 
bulge (217 hours) and 10 for the plane (75 hours), for a total of 292 hours. These 
additional epochs will improve our ability to detect variable sources (but will not 
permit us to conclusively establish the variability phases, at least for many of the 
RR Lyrae stars). These data will also allow the creation of deeper master maps in 
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Table 2: Estimated observing time for the full survey. 



Year 


Time 


RA range 


Moon 


Seeing 


Transparency 


Number of 




[h] 


[hh:mm] 




["] 




K-band epochs 


year 1 


292 


12:00-19:00 


any 


0.8 


clear 


6 (bulge and disk) 


year 2 


292 


12:00-19:00 


any 


0.8 


clear 


4 (bulge and disk) 


year 3 


652 


12:00-19:00 


any 


any 


thin 


80 (bulge only) 


year 4 


525 


12:00-19:00 


any 


any 


thin 


70 (disk only) 


year 5 


168 


12:00-19:00 


any 


0.8 


clear 


20 (bulge) + 9 (disk) 


Total 


1929 













- 6 epochs in Kfor bulge and disk; K |jm =18/20 mag (single/combined epochs) 

- Z,Y,J,H, K single (quasi-simultaneously) epoch observations (bulge & disk) 



- 4 epochs in K for bulge and disk 



ro 3 I 



- main part of bulge variability campaign (80 epochs, 652 h) 

- map bulge and disk once per night 



- main disk variability campaign (similar to bulge, but 70 epochs, 525 h) 



- bulge and disk observations in K band 

- 20/9 epochs spread over the whole year 

- subset will be observed more frequently (1 0-40 times per night) 



t 

Figure 6: Schematic schedule of the VW survey strategy. 



K s , to fine tune the strategy for the main campaign of the following year. 

During the third year, the main bulge-variability campaign of 60 epochs will 
be carried out over 652 hours. According to the ESO Public Survey requirements, 
only 40 nights can be consecutive, while the others will be spread over the bulge 
season. We will use the K s band to map the whole bulge and inner plane. A subset 
of the fields can be observed more frequently (4-8 times per night). This strategy 
allows us to partially remove aliasing and to improve the periods, while being more 
sensitive to shorter timescale variables and microlensing events. 
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During the fourth year, the main Galactic -plane variability campaign will be 
carried out over 70 epochs using 525 hours in K s , following a similar strategy as 
for the bulge in the previous year. According to the ESO requirements, only part 
of these observations can be carried out over 27 consecutive nights. 

Finally, during the fifth year, we will acquire 12 more epochs for the bulge 
and five more for the disk, with observations spread over the season, taking a total 
of 130 and 38 hours, respectively. This allows measurements of longer-timescale 
variables, and the search for high-proper-motion objects. A subset of the fields can 
be observed much more frequently (10^10 times per night). This strategy enables 
detection of short-period variables and planetary transits. 



4 Scientific goals 



The major VVV survey products will be a high-resolution ZY J H K s colour atlas 
of the bulge and plane regions, and a catalogue of variable point sources, includ- 
ing positions, mean magnitudes, and amplitudes. We expect to detect more than 
10 6 variable objects. However, the total number may even reach 10 7 if one takes 
into account results fro m recent deep variable searches in optical bands in se lected 
Galactic-plane regions (iWeldrake & Bayhssll2008uPietrukowicz et all 12009). Our 
database will be public, a significant treasure for the whole community to exploit 
for a variety of scientific programmes. 



The main scientific goals of the VVV survey are: 



1 . To find RR Lyrae in the bulge, which will allow us to determine periods and 
amplitudes, and measure accurate mean K s magnitudes. We will construct the 
Bailey diagram (luminosity amplitude vs. period), and interpret the results of the 
variability analysis in terms of stellar pulsation and evolution mode ls, similarly 
to what is currently done using (primarily) the visual bandpasses (e.g.JBono et al 



200r ll2007l 
2003). As 



2000, 2007; M arconi et al. 



■mg (primarily) the visual bandpasses (e.g., cono et al. 
2003 :lCacciari etall 120051 : iMarconi & DeglTnnocenti 



20071) . As pointed out by iBono et all (12000), one common limitation in the com- 
parison between hydrodynamical pulsation models and the observations is the lack 
of proper sampling in the available near-IR light curves - a situation which will be 
dramatically improved with the advent of the well-sampled A^-band light curves 
from the VVV survey. Naturally, while RR Lyrae stars are the main focus of 
this project, similar such studies will be possible for many other types of variable 
stars, provided they are detected in signi ficant numbers in our studied fields - in- 



cluding, for instance , class ical Cepheids (IFiorentino et all 120071). type II Cepheids 



( Di Criscienzo et al. . 2007 ). anomalous Cepheids ( Fiorentino et al. . 2006T) . and Mi 
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ras and semiregular variables (IKerschbaum et all 120061 : ISmith et all 120061) . among 
others - with a significant impact also upon their use as distance indicators. 

The pulsation properties of bulge variables will be com pared with tho se of 
similar variables in the halo and nearby dwarf galaxies (e.g., ICatelarlE009h . The 
distances measured and RR Lyrae counts can be compared with the red-clump gi- 
ants, which are excellent tracers of the inner bar (iStanek et all 1 19941) . This would 
define the geo metry of the inner bar an d of additional structures (such as a poten- 
tial second bar; Nishivama et al., l2005h . and explore the radial dependence of the 
density (e.g jMinniti et allll999h . or trends with Galactic latitude-longitude, to fi- 
nally unveil the structure of the bulge. The microlensing surveys that we have been 
involved in (OGLE and MACH O) have discovered a small fraction of the bulge 
RR Lyrae stars (e.g., Fig.|7J see Klcock et all 1 19981 : Iwozniak et all l2002h . Our 
survey will increase the amount of data on bulge RR Lyrae significantly. 




2 3 4 

V-l [mag] 

Figure 7: Optical colour-magnitude diagram for fundamental mode RR Lyrae 
stars in the Galactic bulge from the Optical Gravitational Lensing Experiment 
dCollingeetalll2006h 



We will search for RR Lyrae and type II Cepheids in the Galactic bulge to 
reveal the presence of any debris related to accretion events that might have left 
behind the present-day globular cluster NGC 6441. This globular cluster is well 
known to contain an anomalous RR Lyrae population, with periods that are much 
longer than those of known field RR Lyrae stars of similarly high metallicity (e.g., 
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Pritzl et all lioool 2003). In particular, the presence of the unusually long-period 



(P > 0.45 d) RRc (first overtone) variables, which have so far not been found 
in the general fi eld but are present in large numbers in this globular cluster (e.g., 



CatelanL l2004al) . should provide the smoking gun for the presence of NGC 6441- 



related debris in the general bulge field. In a similar vein, long-period RRab stars 
(fundamental-mode pulsators) occupying the appropriate position in the period- 
amplitude diagram should also provide us with a strong indication of prior mem- 
bership of such a protogalactic fragment. 



2. To identify variable stars belonging to known star clusters. There are 33 glob- 
ular clusters and 355 open clusters located in the VVV area (Fig. [8]), which may 
contain RR Lyrae, type I and II Cepheids, semiregular variables, and eclipsing 
binaries, among other types of variables. Distances, reddening values, metallic - 
ities, and horizontal-branch ( HB) types wi l l be obtained for th e se clu sters from 
a homogeneous dataset (e.g., ICatelan et all 120061 : IZoccali et all 120031 ). In some 
favourable cases, ages can be derived. Table[3] lists the globular clusters that will 
be covered, giving positions in equatorial and Galactic coordinates, and distances 
from the Sun. The asterisks in the last column indicate that more than one third 
of these clusters have uncertain distances. We will improve the distances for these 
globulars, and confirm the previous estimates for the rest of the open and globular 
cluster sample. 



-10 -15 Vt -go -it£ l -30 16 -30 15,5 -40 15 n 5-46 H 13 lg.5 IS 




10 -10 -20 -30 -40 -50 

Longitude 



Figure 8: Map of the globular and open cluster positions (full and empty circles, 
respectively) towards the Milky Way bulge and plane. Inc luded in the VVV area 
are the 33 globular clusters listed in Table 3 Jllarris , 1996h . and 355 open clusters 



(|Bica et all 12003k iDutra et all 120031 : iDias et all. 120061) . The bulge contours are 
indicated. 
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Table 3: The 33 known globular clusters within the VVV area. 



RA(J2000) DEC(J2000) I b D 



PliiQtpr IT) 


Th h * m m * c c T 
|_iin. mill. j 


H-IU. . 111111 . f»f» J 






Henri 
|Kpcj 


Terzan 2 


1 T O ^7 O O 1 

17:27:33.1 


-30:48:08 


356.32 


2.30 


8.7 


Terzan 4 


17:30:39.0 


-31:35:44 


356.02 


1.31 


9.1 * 


HP 1 


17:31:05.2 


-29:58:54 


357.42 


2.12 


14.1 * 


Liller 1 


17:33:24.5 


-33:23:20 


354.84 


-0.16 


9.6 * 


NGC 6380 


17:34:28.0 


^ i~\ r\ A A A 

-39:04:09 


350.18 


-3.42 


10.7 


Terzan 1 


17:35:47.2 


o a o n ^ < 

-30:28:54 


357.56 


0.99 


5.6 


Ton 2 


17:36:10.5 


on o o -i o 

-38:33:12 


350.80 


-3.42 


8.1 * 


NGC 6401 


17:38:36.6 


o o ^ ^ o ,( 

-23:54:34 


3.45 


3.98 


10.5 


Pal 6 


17:43:42.2 


-26:13:21 


2.09 


1.78 


5.9 


Djorg 1 


17:47:28.3 


O O AO ^ 

-33:03:56 


356.67 


-2.48 


12.0 * 


Terzan 5 


17:48:04.9 


0,1 /< x A r 

-24:46:45 


3.84 


1.69 


10.3 * 


NGC 6440 


17:48:52.7 


OA O t O ^1 

-20:21:37 


7.73 


3.80 


8.4 


TV T/-1 s~~i f A A t 

NGC 6441 


17:50:12.9 


O i - ? AO A ^ 

-37:03:05 


353.53 


-5.01 


11.7 


TP s 

Terzan 6 


17:50:46.4 


-31:16:31 


358.57 


o 1 y 

-2. 16 


9.5 * 


NGC 6453 


17:50:51.7 


O /( Or" f~7 

-34:35:57 


355.72 


-3.87 


9.6 


UKS 1 


i i—i c A 0^7 O 

17:54:27.2 


O /( AO /IO 

-24:08:43 


5.12 


0.76 


8.3 * 


Terzan 9 


-i O A 1 o o o 

18:01:38.8 


-26:50:23 


3.60 


-1.99 


6.5 * 


Djorg 2 


18:01:49. 1 


0^7 /(A O O 

-27:49:33 


2.76 


O C A 

-2.51 


6.7 * 


Terzan 10 


1 O AO C~l A 

18:02:57.4 


O /" A 1 A A 

-26:04:00 


4.42 


-1.86 


5.7 * 


NGC 6522 


18:03:34.1 


OA AO AO 

-30:02:02 


1.02 


-3.93 


7.8 


NGC 6528 


18:04:49.6 


-30:03:21 


1.14 


-AM 


7.9 


NGC 6540 


18:06:08.6 


-27:45:55 


3.29 


-3.31 


3.7 


NGC 6544 


18:07:20.6 


-24:59:51 


5.84 


-2.20 


2.7 


NGC 6553 


18:09:17.6 


-25:54:31 


5.25 


-3.03 


6.0 


HJK2000-GC02 


18:09:36.5 


_20:46:44 


9.78 


-0.62 


4.0* 


NGC 6558 


18:10:17.6 


-31:45:50 


0.20 


-6.02 


7.4 


Terzan 12 


18:12:15.8 


-22:44:31 


8.36 


-2.10 


4.8* 


NGC 6569 


18:13:38.8 


-31:49:37 


0.48 


-6.68 


10.7 


NGC 6624 


18:23:40.5 


-30:21:40 


2.79 


-7.91 


7.9 


NGC 6626 


18:24:32.9 


-24:52:12 


7.80 


-5.58 


5.6 


NGC 6638 


18:30:56.1 


-25:29:51 


7.90 


-7.15 


9.6 


NGC 6642 


18:31:54.1 


-23:28:31 


9.81 


-6.44 


8.4 


NGC 6656 


18:36:24.2 


-23:54:12 


9.89 


-7.55 


3.2 



uncertain distances 
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For some of the Cepheids which we expect to find in open clusters, we will de- 
rive improved distances by app lying the IR surface-brightness technique (iGieren et al 
2005k iFouque & Gierenl. 1 19971) on these variables. 



3. To find eclipsing binaries in large numbers. We expect to detect more than 
5 x 10 5 binaries, an unprecedented database that will allow us to determine pe- 
riods, amplitudes, mean magnitudes, study stellar properties, and also select ex- 
trasolar planetary-transit candidates. In particular, YY Gem-like systems can be 
identi fied to constrain the lower main-sequence parameters (e.g.. iTorres & Ribas . 
2002T) . and selected transit fi elds can be followed frequently to identify and mea- 
sure extrasolar giant planets (|Udalski et all 12002). 



4. To find rare variable sources. The massive variability dataset and multi-colour 
atlas will allow us to search for CVs (novae, dwarf novae) and other eruptive 
variables (e.g., RS CVn), eclipsing binary RR Lyrae, pre-HB/post He-flash stars 
(ISilva Aguirre et all 120081) . eclipsing binary red-clump giants, luminous blue vari- 
ables (LBVs), FU Ori protostars undergoing unstable accretion, and asymptotic 
giant branch (AG B) stars at the stage of unstable shell burning. The Murchison 
Widefield Array (Lonsdale et all Eo09) will soon open up the field of surveys of 
radio transients in the southern hemisphere, likely making infrared variability sur- 
veys of the southern part of the Galaxy interesting for interpreting their results. The 
Galactic center region contains numerous high energy sources without counterparts 
in the optical/infrared band, with continuing discoveries coming both from space- 
based X-ray and y-ray observatories and from ground-based Cerenkov arrays. The 
most sensitive of these observatories, the Chandra X-ray Observatory, has discov- 
ered almost 10 4 X- ray sources within a 2 x 0.8 deg 2 region near the Ga l actic c enter 
dMuno et all 12009). The Chandra Galactic Bulge Survey dBassa et all 12008) con- 
tains about 10 3 X-ray sources spread over a larger, slightly less crowded region, 
and already has deep optical follow-up. These objects are often extincted enough 
that they cannot be observed in the optical, but there are also often multiple in- 
frared counterparts found even in moder ately deep in f rared images within the 1" 



or smaller Chandra error regions (e.g., iGosling et all 12007). For y-ray sources , 



such as those w hich have bee n found with INTEGRAL (e.g., 



Kuu 



kers et all 2006; 



Bird et all I2007T) and HESS (lAharonian et all 120071 : IChaves et all 120081) . the point 
source localization is even poorer - typically several arcminutes. Fortunately, X- 
ray binaries are usually variable, and t heir variability in i nfrared and X-rays is, 
for the most part, well correlated (e.g., Russell et al. . 2006|) . As a result an in- 
frared variability survey of the inner Galaxy can be combined with the ongoing 
monitoring in X-rays and y-rays from missions like Swift and MAXI, should give 
the opportunity to search narrow down the list of candidate infrared counterparts 
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through searches for IR emitters with variability well correlated with that in the 
high energy bands. 



5. To search for microlensing event, especially reddened events, short time-scale 
events, and high magnification events in obscured high-density fields. The spatial 
dependence of the microlensing opti cal depth t at near-infrared wavelengths has 
been modelled by lKerins et al. (2009) using synthetic population synthesis models 
of the Galactic disk and bulge. The spatial variation of t can probe directly the 
mass distribution contained in the inner regions. Unfortunately, current optical mi- 
crolensing searches do not cover the whole bulge or the plane due to the prevalence 
of dust. In particular, they miss the inner regions of the bulg e where the optical 
depth is higher, providing p oor constraints to current models (IBissantz & Gerhard , 
2002l : lBissantz et alll2004l) . At near-IR wavelengths a map of microlensing optical 
depth for the whole bulge can be made, allowing a search for asymmetries in the 
spatial distribution of t. The strength of the microlensing asymmetry is a function 
of the orientation of the inner bar as well as the relative contributions of the bulge 
and disk to the microlensing rate. A near-IR microlensing sample can therefore 
provide an important additional lever on the 3D geometry of the inner Galaxy. In 
addition, we expect to detect microlensing of stars in the Sgr dSph galaxy (e.g., 
Popowski et alll2005h . 



6. To monitor variability around the Galactic center, an area of 1.5 deg 2 a round 
the Galactic center including the 180 pc Nuclear Ring (IMessineo et all 120021) will 
be the most frequently monitored field, over a total of 200 epochs spanning five 
years. Expected variability due to high- magnification microle nsing, or flares due 
to black-hole accretion, can occur (e.g., IChaname et al. ■ l200lh . Black-hole flares 
easily reach K s ~\6 mag ; with a typi c al du r ation of 10-30 min. T he expected 
flare rate is 2-6 per day (IGenzel et all |2003k Eisenhauer et all 12005), in addition 
to a l onger-timescale variation predicted by accretion simulations (|Cuadra et all 
2006). We also expect some Wolf-Rayet variability in the population of mas- 
sive stars and clusters in this region, and we will search for eclipsing WR stars. 
In addition, we will be able to identify the counterparts of high-energy (y-, and 
X-ray) sources: accreting black holes, microquasars, binary-pulsar companions, 
low mass X-ray binar i es (LMXBs), and high m ass X-ray binaries (HMXBs) (e.g., 
Mirabel & Rodriguez! . ! 19981 : iLucas et al.L l2008h . In particular, this survey may re- 
veal the as yet undetected counterparts of the most luminou s, persistent hard X - 
ray/jet sources in the Galactic center region ; IE 1740.7-2942 ([Mirabel et all 1 19921) 
and GRS 1758-258 dRodriguez et all Il992h . Our survey will allow us to identify 
and monitor the counterparts of several variable high-energy sources, and perhaps 
in some cases to determine their orbital periods. 
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7. To search for new star clusters of different ages and identify their variable 
star members, such as Cepheids, semiregular variables, W UMa-, and 8 Set-type 
stars. The asymmetric distribution of the known globulars in the Gala ctic center re 



gion hints at the presence of additional, as yet undiscovered objects (|lvanov et al 



2005]; iKurtev et all l2008h . Our team members have already carried out success- 



ful campaign s searching for new clusters i n the 2MASS Poi nt Source Catalogue 
dlvanov et all 120021 : iBorissova et"ail 120031 : butra et all l2003h . Note that 2MASS , 
with a Ksj im = 14.5 mag in the bulge, discovered hundreds of open-cluster can- 
didates, as well as two new globular clusters. Because we will reach 3-4 mag 
deeper, we expect to find many new clusters. 

We will build a homogeneous sample of at least 300 open clusters in the di- 
rection of the Galactic center, with accurately derived fundamental parameters. A 
sample thus obtained will be useful to investigate the structure of the Galactic disk 
in directions intercepting the bulge. This will represent a major improvement over 
our current knowledge in this area, since only 20 clusters have so far been studied 
in detail in this region of the sky, according to the last update (ver. 2.10, 2009) of 
Dias' catalogue of Galactic Open Clusters (bias et all l2002h . We will carry out 
a census of the Milky Way open clusters projected onto the central parts of the 
Galaxy. This will allow us to establish the fraction of star clusters compared to 
statistical fluctuations of the dense stellar field in those directions, as well as the 
cluster-formation efficiency relative to field stars. We can also estimate the role of 
disruption effects on timescales < 10 Myr (particularly useful for open clusters). 
We will derive the physical parameters: angular sizes, radial velocities, reddening, 
distances, masses, and ages of these clusters. Reliable fundamental parameters 
of unstudied open clusters are important both to disk studies and to constrain the 
theories of molecular-cloud fragmentation, star formation, as well as stellar and dy- 
namical evolution. We will trace the structure of the Galactic disk. Recent studies 
of the disk structure based on open clusters are complete up to only 1 kpc from the 
Sun. We will complement and re-derive the existing kinematics distributions such 
as distance of the cluster to the Galactic center vs. age distribution, open clusters 
age histogram, distance of the open clusters to the Sun vs. reddening, etc. 



8. To provide complementary near-IR multi-colour information (reddening, tem- 
peratures, luminosities) and time coverage to the following past and on-going sur- 
veys: GLIMPSE-II, VPHAS+, MACHO, OGLE, EROS, MOA, Pan-STARRSl, 
and PLANET. Near-IR photometry is important for the events discovered by mi- 
crolensing surveys. For old previously detected events or new ones, the VVV 
survey will give field reddening, a baseline colour, and a magnitude that can im- 
mediately be translated into temperature and luminosity for the source star. Char- 
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acterization of the source is essential for refini ng the microlensing l ight-curve pa- 
rameters and the physical lens properties (e.g. jBeaulieu et al.L I2006Q . 



9. To find variable stars in the Sgr dSph galaxy: Figure [7] shows that the Sgr dSPh 
RR Lyrae are well within reach and can be r eadily i dentified. RR Lyrae woul d 
give the 3-D structure of the Sgr dSph (e.g., klardl 1 19961 : Klcock et all Il998u . 
To measure the depth and the tilt of the Sgr dSph along the line of sight, mean 
RR Lyrae magnitudes with an accuracy better than 0.01 mag are necessary. This 
corresponds to the t idal radius of M54, supposedly the core of the Sgr dSph galaxy 
(IMonaco et all 120041) . We will also detect and measure carbon stars, semiregulars, 
and eclipsing binary members of the Sgr dwarf galaxy (type II Cepheids are ex- 
pected in the bulge, but classical Cepheids are not; since these would be found in 
the disk instead). 



10. To identify high-proper motion objects and background Quasi-Stellar Objects 
( QSOs): this goal links the - seemingly unrelated - intrinsically faintest and bright- 
est objects in the Universe. On the faint end, we would use proper motions to find 
nearby late M-type stars, brown dwarfs (L and T types), and high-velocity halo 
stars. The proper motions of the most interesting low-mass objects will be deter- 
mined using UKIDSS, DENIS, and 2MASS to extend the time baseline, and to 
search for objects with smaller proper motions. With an astrometric accuracy o f 
~ 10 mas, comparable to UKIDSS (e.g- Lodieu et all 120071 : beacon et al.L l2009h . 
or slightly lower in the most crowd ed regions, we will be able to determine this 
proper motions ( Lucas et al. . 20081) . Scaling the results of the UKIDSS proper 
motion studies we expect an accuracy of < 15 miliarcs yr _1 . This is beyond the 
minimum value to distinguish bulge and disk populations, since the proper mo- 
tion difference between them is ~ 6 miliarcs yr~' in the galactic center region 



(IClarkson et al.ll2008l) . 

On the intrinsically bright end, variability would also allow us to identify back- 
ground qu asars, provid i ng an extragalactic reference scale for future proper mo- 
tions (e.g., IPiatek et all 120051) . QSOs have a relatively broad colour range depend- 
ing on their redshif t, and their intrinsic v ariability increases monotonically with in- 
creasing tim e lags ( de Vries et all 120050 . Their amplitudes should be > 0.2 mag in 
the near-IR (lEnya et all 120021) . We estimate that we will find > 500 active galac- 
tic nuclei (AGNs), assuming a surface density of 2 deg with K s < 15.5 mag 
(ILeipski et all 120051) in the regions above and below the disk where Ak < 0.5 mag. 



11. To identify pre-main-sequence (pre-MS) clusters and associations through 
variability: neither the duration of pre-MS evolution as a function of mass nor 
the duration of active star formation in molecular clouds are well established. IR 
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studies (e.g. ICarpenter et all 120011) have shown that a large fraction of all pre-MS 
stars are variable over a wide range in absolute magnitude, therefore, the identifi- 
cation of young stellar objects (YSOs) via a variable survey is a suitable method 
to detect YSOs due to its independence on color. VVV will be capable of picking 
out loose associations of pre-MS stars in crowded fields long after the molecular 
cloud has dispersed and the cluster has become unbound. Through careful anal- 
ysis of the detection rate of such dispersed pre-MS populations and comparison 
with the detection rate of pre-MS accretion disks with Spitzer and WISE it will be 
possible to determine the relative durations of the phases with and without a disk. 
Two such rare pre-MS variables, KH 15D and VI 648 Ori, have been de scribed by 



Herbst et all (|2002h : lKusakabe et all (|2005D and lReipurth & Aspinl (120041. 



5 Conclusions 

We have described our near-IR public survey of the inner Milky Way bulge and 
disk. The VVV is a survey to be carried out with VISTA at Paranal Observatory 
between 2010 and 2014. It will map repeatedly most of the Milky Way bulge, as 
well as the inner southern disk, covering a total area of about 520 deg 2 containing 
about ~ 10 9 point sources, 33 known globular clusters and more than three hun- 
dred known open clusters. The main survey products will be a Z F J H K s atlas of 
the Milky Way bulge and inner disk, and catalogues of variable point sources and 
high-proper-motion objects. The multi-epoch i^-band photometry will allow the 
identification and phasing of periodic variable stars, as well as microlensing events 
and planetary transits. We will unveil the 3-D structure of the inner bulge and disk 
of the Milky Way using well understood distance indicators such as RR Lyrae stars, 
Cepheids, and red-clump giants. The survey is expected to detect more than 100 
star forming regions and pre-MS clusters, benefiting from sensibility to variable 
pre-MS stars in regions where the molecular cloud has been dispersed. This large 
sample will permit measurement of the duration of the pre-MS evolution as a fuc- 
tion of mass. It will also allow investigation of the effect of environment on the 
outcome of the star formation process. The VISTA observations will be combined 
with data from MACHO, OGLE, EROS, 2MASS, DENIS, HST, Spitzer, Chandra, 
INTEGRAL, AKARI, WISE, Fermi LAT, XMM-Newton and in the future GAIA 
and ALMA for a complete understanding of the variable star sources in the inner 
Milky Way. 
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